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1) ' 
^G ' Abstract. We introduce a natural model of quintessence in string theory where the light 

rolling scalar is radiatively stable and couples to Standard Model matter with weaker-than- 

Planckian strength. The model is embedded in an anisotropic type IIB compactification with 

i^T two exponentially large extra dimensions and TeV-scale gravity. The bulk turns out to be 

\/~} • nearly supersymmetric since the scale of the gravitino mass is of the order of the observed 

^^O , value of the cosmological constant. The quintessence field is a modulus parameterising the 

size of an internal four-cycle which naturally develops a potential of the order (gravitino 

£Lj ' mass) 4 , leading to a small dark energy scale without tunings. The mass of the quintessence 

{vq ■ field is also radiatively stable since it is protected by supersymmetry in the bulk. Moreover, 

this light scalar couples to ordinary matter via its mixing with the volume mode. Due to 

the fact that the quintessence field is a flat direction at leading order, this mixing is very 

small, resulting in a suppressed coupling to Standard Model particles which avoids stringent 

fifth- force constraints. On the other hand, if dark matter is realised in terms of Kaluza-Klein 



states, unsuppressed couplings between dark energy and dark matter can emerge, leading to a 
scenario of coupled quintessence within string theory. We study the dynamics of quintessence 
in our set-up, showing that its main features make it compatible with observations. 
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1 Introduction 

Present cosmological observations indicate that our Universe is currently accelerating. This 
accelerated expansion is driven by a mysterious form of energy density with negative pressure, 
dubbed 'dark energy', that accounts for about 70% of the total energy density of the Universe. 
The major part of the remaining 30% consists of cold non-baryonic dark matter, whose 
microscopic origin is not known, although it is likely to be associated with new particles as 
supersymmetric partners, axions, string moduli or hidden sector degrees of freedom. Finally, 
ordinary baryonic matter accounts for only a few% of the total energy density. 

The microscopic origin of dark energy is even less clear than that of dark matter. Many 
proposals have been pushed forward for explaining present day acceleration (see [1] for a 
comprehensive review). We will focus here on quintessence models, in which the dynamics of 
scalar fields drives the present cosmological acceleration, as opposed to the simplest scenario 
in which the cosmological energy density is dominated today by a pure cosmological constant. 

Both of these proposals, a rolling quintessence field and a non-zero cosmological con- 
stant, have their own virtues and shortcomings: 

• Cosmological Constant Problem: If dark energy is given by a non-zero cosmological 
constant, in order to match the current cosmological observations, its value needs to 
be of the order (V) = A 4 ~ meV ~ 10 -120 M p . Whereas, if dark energy is provided 
by a quintessence model, in which a late-time acceleration is characterised by a scalar 
(f) with kinetic energy smaller than its potential energy, the corresponding scalar mass 



has to be of the order m ~ 10 33 eV. The two values of A and m are related to each 
other since: 

™ 2 ~ V H ~ ^ Jj2 ~ 10" 120 M| ~ (l(T 33 eV) 2 . (1.1) 

The condition 4> ^ Mp follows from the Friedmann equation 3MpH 2 = ((p 2 /2+V) ~ V 
and the equation of motion for (j), = (f> + 3H(p + V& — 3i70 + 16,. In fact, for (p 2 < V, 
we have: 

^ ~ T7" ~ T7T £ M ^ ■ ( L2 ) 

The required smallness of these quantities (A or m, depending on the scenario under 
consideration) represents one of the biggest mysteries of modern physics since a naive 
estimate of the vacuum energy would yield a result many orders of magnitude larger 
than the observed one: (V) ~ M p for non-supersymmetric theories and (V) ~ TeV ~ 
10 -60 M p for theories with TeV-scale super symmetry. In the quintessence case, one 
can assume that the vacuum energy is exactly zero (as forced by symmetry arguments, 
for example), but one should still explain how to obtain such a small mass for the 
rolling scalar. Indeed, scalar masses are notoriously hard to keep from getting large 
contributions when integrating out ultra-violet physics running in the loops. This is 
the issue of technical naturalness, which plays a crucial role in the determination of the 
Higgs mass and the solution of the famous ^-problem in slow-roll inflation. However, 
here the problem is more severe than in inflation since one has to integrate out all 
the modes down to the tiny cosmological constant scale A ~ O(meV), instead of the 
inflationary scale M; n f ~ O(Mqut)- 

Coincidence Problem: Observations indicate that the energy densities of dark matter 
and dark energy are of the same order of magnitude at present cosmological epochs. 
The reason why different forms of matter have similar properties today is a puzzle which 
goes under the name of 'coincidence problem'. Dark energy models based on a non-zero, 
pure cosmological constant have little to say about it, but quintessence models can in 
principle provide answers to this problem. The idea is to consider 'tracking' solutions, 
in which the non-trivial dynamics of the scalar field leads its energy density to mimic 
radiation during the radiation-dominated era. The transition between radiation to 
matter domination, then, triggers quintessence towards assuming characteristics that 
match the observed behavior of dark energy. 

Fifth-force constraints: Despite providing a potential understanding of the present 
ratios of energy densities via tracking solutions, quintessence models are generically 
plagued by the problem of long-range unobserved fifth-forces which would be mediated 
by such a light scalar field [2] . The present mass bounds for scalars with Planck strength 
couplings to matter are m > O(meV) [3]. Furthermore, a rolling scalar would lead to 
a time variation of the constants of Nature, while these do not appear to drastically 
change between Big Bang Nucleosynthesis (BBN) and the present epoch. The standard 
solution to this problem involves pseudo Goldstone bosons with shift symmetries [4] 
since these pseudoscalars give rise to spin-dependent couplings which cannot be probed 
in fifth-force experiments that are insensitive to the polarization. Moreover, the shift 
symmetry forbids higher dimensional operators, rendering the quintessence field stable 
against radiative corrections from UV physics. On the other hand, it is much harder 
to satisfy the current constraints in models where the quintessence field is a scalar. 



Most of the attempts to realise quintessence within string theory rely on axions and shift 
symmetries [5-8]. In this paper we shall, however, take a different approach, and build the 
first string model where the quintessence field is a scalar parameterising the size of an internal 
four-cycle. We shall focus on LARGE Volume Scenarios (LVS) which emerge naturally in 
type IIB flux compactifications. The set-up of this model represents an interesting variation 
of the one derived in [9] for inflationary purposes. In fact, as in [9], we focus on the recent 
string embedding [10] of the Supersymmetric Large Extra Dimension (SLED) proposal [11]. 
The field which was playing the role of the inflaton in [9] will now behaves as quintessence. 

In our new model, quintessence is driven by a closed string modulus with a naturally 
small mass m protected against loop corrections. The mechanism for obtaining such a light 
scalar whose mass is radiative stable is similar to the one of [12] and we shall outline it in 
section 2. Contrary to the model of [12] where the quintessence field was the overall volume 
mode, in our case quintessence is driven by a fibre modulus that at leading order neither 
develops a potential nor couples to Standard Model (SM) matter. The potential for the 
quintessence field arises only at subleading order by means of non-perturbative and string 
loop effects, that also induce suppressed, weaker-than-Planckian couplings with ordinary 
matter. Thus, our model is able to avoid fifth-force constraints. An intuitive explanation for 
this non-standard behaviour for the quintessence scalar will be given in section 2, while in 
sections 3 and 4 we will develop our arguments at a more technical level. 

The dynamics of the quintessence-dark matter system will then be analysed in section 5 
starting from the standard decoupled case and then moving to the more complicated coupled 
case where we shall also explain how our stringy set-up can allow for unsuppressed couplings 
of the quintessence field with dark matter degrees of freedom. Hence we will suggest a 
realisation of a coupled quintessence-dark matter model inspired by string theory. Finally 
we will conclude in section 6. 

2 Natural quintessence from anisotropic compactifications 

In this section, we shall outline in some detail the ideas behind the construction of our 
quintessence model, before developing the technical elements of our set-up in the next sec- 
tions. 

2.1 Naturally light quintessence field 

In our approach, a naturally small mass for the quintessence field is generated by means of a 
string realisation of the so-called Supersymmetric Large Extra Dimension (SLED) scenario 
(see [13] for a review). In its original form, this is a brane- world model with two large 
extra dimensions, in which supersymmetry is preserved in the bulk, whereas it is only non- 
linearly realised on a Standard Model codimension-two brane. Its main attractive feature 
is the possibility to tie the discussion of the hierarchy problem to that of the cosmological 
constant problem. More precisely, the hierarchy problem is addressed by the presence of 
large extra dimensions which lower the fundamental gravity scale down to the TeV-scale; 
at the same time, the observed value of the cosmological constant is associated with the 
low supersymmetry breaking scale in the bulk. These two arguments are tied up since 
simple considerations based on dimensional reduction yield the following relation among the 
fundamental (4 + (i)-dimensional gravity scale M* and the gravitino mass m 3 / 2 : 

™3/2~h^J M P . (2.1) 



Setting M* ~ O(TeV) ~ 10~ 15 Mp, the gravitino mass turns out to be very small, of the 
order of the dark energy scale m, 3 / 2 ~ A ' ~ O(meV) ~ 10~ 30 Mp. This implies that the 
bulk is nearly super symmetric, given that supersymmetry is broken only at the dark energy 
scale. Quantum corrections to the vacuum energy from loops of bulk fields scale as: 

Moo P ~ Str(M 2 ) A 2 VV ~ m 2 /2 M 2 K , (2.2) 

where Str(M 2 ) ~ m \/2 an( ^ * ne cut-off is taken to be the scale at which a pure 4D description 
ceases to be valid: Auv ~ Mkk- Given that the Kaluza-Klein scale is given by: 

M KK ~ l-j^-J m p . ( 2 - 3 ) 

this is of the same order of m 3 / 2 for d = 2. This is also the case that for TeV-scale fundamental 
gravity leads to micron-sized extra dimensions which are at the edge of detectability in 
experiments that look for deviations from Newton's law. Hence in the case of two large extra 
dimensions, the loop generated potential (2.2) gives the correct observed order of magnitude 
of the cosmological constant since it scales as V\ oop ~ m 3/2 ~ 10 -120 Mp. Notice that for 
energies above the cut-off Mkk ; the theory is effectively 6D and supersymmetric, and so the 
contributions to the vacuum energy from loops of bosons and fermions cancel among each 
other. 

In order to have a viable and complete solution of the cosmological constant problem, 
one has to make sure that (2.2) is the leading order contribution to the vacuum energy. 
In general, one would expect larger contributions coming directly from the tension of the 
Standard Model brane, of the order T ~ O(TeV) 4 , and from the tree-level part of the bulk 
potential corresponding to the curvature of the extra dimensions induced by the presence of 
the TeV brane. This is a delicate issue (see e.g. [14]): recent progresses seem to suggest that 
consistent scenarios can be built along these lines by means of suitable couplings between 
bulk fields and the SM brane [15]. Since the rest of our discussion does not rely specifically 
on these details, we will not develop this issue any further in this paper. 

The SLED framework opens the possibility to obtain a natural model of quintessence. 
In fact, in 6D SLED models, the radion mode r develops a potential via loops of bulk fields 
(Casimir energy) which scale as (2.2), and so give it a mass of the order [12]: 



m„ ~ 



2 Hoop m 3/2 hn-33 \A 2 Ic, A \ 



Mp M p 

This is the perfect scale for quintessence, since it is the typical mass scale expected for a 
quintessential scalar; for the arguments we outlined above, this small mass is technically 
natural since it is radiatively stable thanks to the unbroken supersymmetry in the bulk (see 
[12] for more details). 

In the quintessence model based on the 6D SLED scenario, the radion mode, being a 
scalar instead of a pseudo-scalar with Planck strength coupling to ordinary matter, would 
mediate long range fifth- forces. Moreover its cosmological evolution would lead to a time 
dependence of Newtons constant. The authors of [12] propose to cure this problem exploiting 
the fact that the couplings of the radion to matter are field dependent, approaching small 
values for appropriate cosmological evolutions of the radion field. The strongest constraint 
comes from the requirement that the Newton's constant does not change that much from 
BBN till today. 



We will now argue that a natural, automatic solution to the fifth-force problem emerges 
when embedding this quintessence scenario in string theory. A recent paper [10] discusses 
how to embed these 6D SLED scenarios within type IIB flux compactifications. The simplest 
Calabi-Yau three-fold that allows such a derivation has a fibred structure with a 4D K3 or 
T 4 fibre over a 2D IP 1 base. This structure allows to take a very anisotropic limit where the 
base is exponentially larger than fibre, resulting in an intermediate 6D effective field theory 
when compactifying 10D type IIB string theory down to 4D. The authors of [10] show how 
to stabilise the geometric moduli in such a way to obtain the above mentioned anisotropic 
limit within the framework of the LARGE Volume Scenario (LVS) [17]. In this framework, 
exponentially large extra dimensions emerge from a moduli stabilisation mechanism that 
exploits a combination of non-perturbative and ol effects. This set-up provides the perfect 
arena for building a SLED-motivated quintessential scenario in IIB string theory: in the next 
subsection we will outline why and how it automatically solves the problem of fifth- forces. 

2.2 Naturally decoupled quintessence field 

The explicit string embedding reveals that the 6D effective field theory is far more richer 
than the one considered in [12], due to the presence of additional closed string moduli on 
top of the radion (that in the IIB embedding corresponds to the volume mode). Moreover, 
besides the equivalent of the loop corrections (2.2), both higher-derivative contributions (c/ 
corrections) and non-perturbative effects play a crucial role in determining the bulk action. 
These additional features of the 6D effective field theory, obtained as the low-energy limit 
of an anisotropic fibred compactification, might at first glance seem a source of new compli- 
cations for getting a viable quintessence model. However, a deeper analysis reveals that the 
string embedding actually provides an elegant quintessence model, with a natural solution 
to the fifth-force problem. 

The minimal scenario we consider in this work involves three Kahler moduli, which 
parameterise the size in string units of internal four-cycles, and a fibered Calabi-Yau volume 
of the form [16]: 

V = t b T f -T^ 2 . (2.5) 

The fields defining the volume have the following roles and geometrical interpretation: 

• t s is a del Pezzo divisor supporting non-perturbative effects of the form W np = A e~ aTa , 
where r s = Re(T s ) and A and a are r s -independent constants. These effects fix r s at 
'small' values but still in a geometric regime above the string scale: (r s ) ~ O(10). 

• The overall volume mode V — £&7y is a modulus fixed at exponentially large values, 
(V) ~ e 1 ' 9 * ^> 1 for weak string coupling g s < 0.1, by the interplay of non-perturbative 
effects and a' corrections to the Kahler potential [17]. The size of the 2D P 1 base is 
given by t b whereas the volume of the 4D K3 or T 4 fibre is given by tj. The value of 
V plays a crucial role in determining the scales and coupling strengths of our model, 
explaining in a natural way why the size of dark energy is so small. 

• The fibre modulus tj plays the role of the quintessence field. It develops a potential only 
a subleading order due to poly-instanton corrections to the superpotential generated 

/ , where 
again rt = Re(Ty) and C is a -ry-independent constant. These tiny effects fix the fibre 
modulus 'small', i.e. at the same order of magnitude of t s : (t/) ~ (r s ). This gives rise 
to a very anisotropic compactification with £& exponentially larger than -JrJ. 



None of these three four-cycles can support the visible sector where supersymmetry is non- 
linearly realised via a proper configuration of branes wrapped around internal divisors. 1 
In fact, the Standard Model cannot be localised neither on t s nor on tj due to the tension 
between chirality and the generation of non-perturbative effects [18]. Moreover, the four-cycle 
containing the base, r& ~ t^y/fj, would give rise to an exponentially small gauge coupling 
g~ 2 ~ r& S> 1. Thus we need to add a fourth Kahler modulus for the Standard Model, tsm- 
This divisor has to be rigid in order to avoid undesired matter in the adjoint representation 
and has to be fixed at small values in order to get a correct value of the gauge coupling. 
This can be done by exploiting D-terms and string loop corrections to the Kahler potential 
[19]. Therefore tsm has the same features as r s apart from the fact that one cycle is fixed 
perturbatively whereas the other non-perturbatively. Hence, without loss of generality, we 
shall just restrict ourselves to the simplest case with three moduli having in mind that the 
modulus t s can represent a generic rigid and 'small' divisor which can be either tsm or the 
one supporting the non-perturbative effects. 

The structure of the scalar potential for the bulk fields (trading £& for V) is the following: 

V = V np (V, Ts) + Vaf (V, Ts) + V poly (V, T s ,Tf) + Moop (V, T s ,T f ) , (2.6) 

where the leading order effects are given by the non-perturbative and a' potentials, V np and 
V a i, which depend only on V and t s and scale as: 

Mead = V np (V, T s ) + Vd (V, T S ) ~ O (rojj /2 M P ) . (2.7) 

The previous potential Vi ea d determines the properties of the volume and the modulus r s at 
leading order in an inverse volume expansion (since the volume is exponentially large, its 
inverse represents an ideal expansion parameter). This implies that the basic features of SM 
physics, like gauge and Yukawa couplings, are determined at leading order by the VEV of 
T a - Instead, in our string construction, the fibre modulus Tf develops a potential only at 
subleading order due to poly-instanton effects and loops of closed strings propagating in the 
bulk. The subleading potential scales as (2.2) : 2 

V suh = V poly (V, T a ,T f ) + V loop (V, T s ,y) ~ O (mf /2 ) . (2.8) 

If the tension of the Standard Model brane is compensated by the brane back-reaction, as in 
the 6D SLED scenario, the leading order contribution to the vacuum energy is given by (2.7) 
which does not reproduce the correct value of the cosmological constant for m 3 / 2 ~ O(meV) 
obtained from (2.1) setting the string scale M s = M* of the order the TeV scale. However the 
potential (2.7) gives rise to a SUSY breaking AdS vacuum, which still needs to be uplifted to a 
dS solution. This can be successfully accomplished by means of one of the many mechanisms 
known in the literature (tension of anti branes at the tip of warped throats [20], D-terms 
from magnetised branes [21] or dilaton-dependent non-pcrturbativc effects [22]). The leading 
order contribution to the vacuum energy would then be given by the subleading potential 



1 We need to focus on a non- linear realisation of supersymmetry where the low-energy effective field theory 
contains just the Standard Model without any superpartner since models with TeV-scale strings predict m 3 / 2 ~ 
meV, and so gravity mediation would give rise to soft terms of this unacceptably small order of magnitude. 

2 We do not include the contribution coming from loops of open string fields since by construction no 
D7-brane stack wraps either the fibre or the base divisor. Hence open strings are localised far away in the 
Calabi-Yau, and so do not develop any potential for n. If that were the case, the internal manifold would 
have a isotropic shape since r\ would also be fixed at exponentially large values: (n) — gj (V) 2 ' 3 S> 1. 



(2.8) which now naturally gives rise to the observed size of dark energy, since it scales as the 
bulk loop potential of 6D SLED models (see (2.2)). 

A key difference between our scenarios and SLED constructions is the fact that the 
volume mode can no longer be used as the quintessence field, since it develops a potential at 
order m%, 2 Mp, and so the corresponding mode is too heavy to drive quintessence: 

2 Mead 2 f m 3/2\ / ln -18 ,rt 2 fo n\ 

m v ~ -^T ~ ™ 3 /2 \j£ ) ~ (10 eV) . (2.9) 

On the other hand, the fibre modulus r/ does not appear in the expression (2.7) for VJ ea d, and 
so its mass has the perfect quintessential value which can be estimated from the subleading 
potential (2.8): 

^~W~ ro ^) 2 ~ (10 " eV)2 ' (2 ' 10) 



This expression is similar to (2.4) which gives the mass of the radion mode of vanilla 6D 
SLED models. Consequently the natural explanation of why the quintessence field is light 
and radiatively stable is, in practice, the same in both scenarios. However, the radion mode 
in the SLED proposal couples to ordinary matter with Planck strength, 1/Mp. Instead, 
as discussed also in [10], the fibre modulus has a weaker-than-Planckian coupling of the 
order e/Mp where e ~ m^/ 2 /Mp <C 1. For this reason, this set-up is ideal for building 
a quintessence model since the rolling field can safely evade all the current bounds from 
fifth- force experiments. 

However, given that V couples to ordinary matter as 1/Mp, these bounds would not 
be evaded by the volume mode if it had a mass below the meV scale as estimated in (2.9). 
Nonetheless, supersymmetry is only non-linearly realised on the SM brane, and so the moduli 
masses (2.9) and (2.10) might get large radiative corrections from loops of open strings. 3 As 
estimated in [10] , this is indeed the case for the volume mode whose mass receives a correction 

of the order: 

M 2 
Smv-j^-- m 3/2 ~ O (meV) , (2.11) 

which renders this modulus still compatible with fifth-force bounds. We stress that the mass 
of the quintessence field does not get lifted by loops of brane fields due to the weakness of 
its coupling. In fact, the correction estimated in [10] scales as: 

M 2 

8m Tf ~ e^- ~ em 3/2 ~ O (HT 33 eV) , (2.12) 

and so does not lift the mass (2.10). 

The intuitive reason of the weakness of the fibre modulus coupling to ordinary matter 
is the fact that the leading order potential Vjead does not depend on tj (see (2.7)), and so 
there is no leading order mixing between ry and the volume mode V. The mixing between 
Tf and V arises only at subleading order at the level of V su \> (see (2.8)). Therefore the order 
of magnitude of this small mixing, e ^C 1, is weighted by the ratio between V su h and V\ ca d- 

Vlead V M P J 



3 These loops of brane fields contribute to the brane tension which has to be cancelled by the brane back- 
reaction. However this cancellation concerns just the vacuum energy, i.e. the value of the potential at the 



minimum. 



Given that the volume mode couples to ordinary matter as 1/Mp and the coupling of the 
fibre modulus to SM particles is induced by the mixing between Tf and V, the coupling of 
the fibre modulus has to be suppressed with respect to the one of the volume mode by the 
small mixing parameter e, explaining the weaker-than-Planckian nature of this interaction. 
In other words, the weakness of this coupling is reflecting the fact that Tf is a flat direction at 
leading order and does not couple with ordinary matter at that level. On the other hand, the 
quintessence field has the opportunity to couple with Planckian strength to massive degrees 
of freedom (as stable 6D Kaluza-Klein states) that can provide a model of dark matter. 
Consequently, our set-up can lead to a stringy realisation of an interacting dark energy-dark 
matter system. 

We shall make this analysis more precise in section 4 where we will compute the moduli 
couplings to Standard Model particles. These are extracted from the moduli dependence of 
the kinetic and mass terms of gauge bosons and fermions, using the canonical normalisation 
carried out in appendix A. 

After this outline of the properties of our quintessence scenario, we pass to develop the 
technical details. 

3 Explicit form of the quintessence potential 

In this section, we start developing our set-up by discussing the features of the potential that 
drives quintessence. As explained in section 2, we work within the framework of the LARGE 
Volume Scenario of type IIB string theory [17] where the volume is fixed at exponentially 
large values by combining a' corrections to the Kahler potential K and non-perturbative 
contributions to the superpotential W. 

As explained in subsection 2.2, we consider compactifications where the volume of the 
Calabi-Yau manifold has a fibered structure and is given by (2.5). With this choice of internal 
manifold, the standard LVS framework stabilises only two of the three Kahler moduli. The 
remaining flat direction, which can be identified with the fibre modulus Tf, is lifted either 
by string loop corrections to K [23] or by poly-instanton effects [9, 10]. Here we focus our 
attention on the latter possibility, and investigate if the resulting potential for the fibre 
modulus Tf allows us to provide a string theory realisation of quintessence. The analysis 
closely follows that of our previous paper [9], in which we studied a similar scalar potential 
in the context of an inflationary model. For this reason, we will only sketch the derivation 
and the properties of the scalar potential here, referring the reader to [9] for more details. 

It was shown in [24] that there are instantons which do not give rise to a single contribu- 
tion to the superpotential but, due to the presence of extra fermionic zero-modes, they correct 
the action of another instanton wrapping a different internal cycle. These poly-instantons 

generate corrections to the superpotential of the form W np ~ e~\ i+e > . These double 
exponential terms were used in [10] to provide a string realisation of 6D SLED scenarios and, 
in [9], to provide a controlled model of inflation in the Kahler moduli sector. 

In what follows we assume the following brane set-up: the cycle t s is wrapped by 
a stack of magnetised D7-branes where appropriate gauge fluxes give rise to two different 
gauge groups that separately undergo gaugino condensation. This gives rise to a race-track 
superpotential. Due to the existence of a Euclidean D3- instanton wrapping the fibre Tf, 
this racetrack superpotential acquires further corrections. Following [24], the gauge kinetic 
functions of the two condensing gauge groups receive non-perturbative corrections, resulting 



in a poly-instanton corrected superpotential of the form: 

W = W + ie- a ( T ' +Cie " M/ ) - Be<™*~™') , (3.1) 



Here, A and B are threshold effects, while a = 2ir/n a and b = 271/11^ with n a ^ natural 
numbers. C\ and C2 are constants associated with the effects of the D3-instanton wrapping 
Tf. On the other hand, the 0/ corrected Kahler potential reads: 

K = K + K a , ~-21nV--g|-, (3.2) 

9s V 

where £ is a positive topological quantity. 

The scalar potential computed from eqs. (3.1) and (3.2) can be written as: 

V = Vb (v -3) (r s , V) + V 0(v -z- P) (r f ,T s , V) , (3.3) 

where the first term scales with the volume as V~ 3 , while the other as V~ 3 ~ p , with p a 
positive, order one parameter. If V is large, the second piece is much suppressed with respect 
to the first one: the previous equation reproduces the split between leading and subleading 
potentials we discussed in section 2 (see eq. (2.6)). The leading contribution to the scalar 
potential, V@(y-3\, is the standard LVS potential that stabilises the blow-up modulus r s and 
the volume mode V: 



V, 



0(V" 3 ) 



8^ (A 2 a 2 e- 2aT ° + B 2 b 2 e~ 2bT ° - 2AB abe 

3V 
4W t s (Aae~ aTs - Bbe~ bT °) 3£W 2 



-(a+b)r 3 



(3.4) 



V 2 A*/ 2 



4<£' V 3 



The subleading Vo(y-3- P \ term depends on all three directions of Kahler moduli space and 
generates a minimum for the fibre modulus Tf. 



16Ci^e- 2brs 



V 0{v - S - P) = - e- 2 ^f \ - 1DOlV o ;f [Z 2 a - ZB(a - b)] 



-br 3 



r 4Cl ^f - [2tt Zr f + Zar s - Bb(a - b)r s } (3-5) 



V 2 

+n 



16Bb 2 '\pT s e- 2hTa Z 4WnB&e -&r ' 



+ "'\n -(br s + 27TT f ) 



3V V 2 

where Z = Bb - Aae'^-^ and n = C 2 -C 1 . 

Given that the mass spectrum in the Kahler moduli sector exhibits the following hier- 
archy m 2 ~ V~ 2 ~> m 2 ) ~ V -3 » m 2 r ~ V~ 3_p with p ~ 0(1) [9], we set r s and V to their 
minima determined (analytically) by minimising the potential (3.4): 

{rs?' 2 = ^ m , (V) = h e b ^ ■ (3.6) 

32a 7 / 1 (l-/ 1 ) 5s 3/2 

Here we have defined: 

(a-b)r s Bb + Z(l-aT s ) 3 ai W ^j 

h = 77 m — uh 1 7/1 5 V ' and h = 7 /1 • ^ 3 - 7 ) 

4 (a — 0) r s i>o + z(l — 4ar s J z 



This allows us to reduce the analysis of quintessence to the study of a single field potential 
for Tf, the subleading potential of eq. (3.5) where we substitute the VEVs of t s and V. 
Following this procedure and performing the canonical normalisation of the fibre modulus 
(see appendix A, in particular the discussion around eq. (A. 9)), the poly-instanton generated 
potential is written as: 



1/ -^poiy 

K P ol y (v) 3 +p 



;| e s/3 _1 ) / / 20 



1 + c eVa 



(3.8) 



where <j> is the field parameterising the shift of the canonically normalised fibre modulus from 
its minimum: <j) = <fi — (eft) = \/3/21n(l + fj/(ry)). The constant c is given by c = 27t(tj) = 
pb{r s ) + 1. K>oiy, on the other hand, is given in terms of the underlying parameters as: 

F poly = fa W »f^ h T\e-iW^Z-nBb). (3.9) 

Equation (3.8) provides a first, extremely flat contribution to the quintessence potential we 
will be working with. Note that, as explained in section 2, we included in eq. (3.8) an 
uplifting term, V^ p = F po i y /(V} 3+p so to set the minimum of the quintessence potential to 
zero. 

As pointed out in [9], there will be further contributions to the fibre modulus potential, 
in particular UV physics will add new <j> dependent terms to the total potential. These are 
loop contributions that we have to take into account, since they can influence the dynamics 
of the quintessence field. These contributions can be estimated via the Coleman- Weinberg 
potential and are due to subleading g s corrections to the Kahler potential. Supersymmetry 
reduces their size, and following the straightforward estimate performed in [9], we obtain: 

Vioop - ^ e vV , where F loop = |- (g s C loop W ) 2 . (3.10) 

The full potential for the fibre modulus that corresponds to our quintessence field is then 
given by: 

Vfibre = Vpoly + Moop, (3-H) 

which we plot in figure 1. Given that the gravitino mass in terms of the overall volume is 
m 3 / 2 ~ Mp/V, this potential for p ~ 0(1) scales as Vabrc ~ (Mp/V) ~ rn t/2^ reproducing 
the behaviour described in section 2. 

For the study of the dynamics of the quintessence field it is convenient to divide the 
potential into three distinct regions: 

1. The region I corresponds to the vicinity of the minimum where the potential is very 
steep; 

2. The region II is the poly-instanton dominated part of the potential; 

3. In region III the loop corrections are the dominant contribution to Vabrc- 

Notice that, as outlined in section 2, we have obtained an extremely flat quintessential poten- 
tial (above all in region II) in which supersymmetry in the extra-dimensions allows to tame 
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Figure 1. Quintessence potential for the fibre modulus. The dotted line represents the poly-instanton 
generated potential, the dashed line the loop potential from UV physics and the full line the total 
potential for the canonically normalised fibre modulus. 



the dangerous corrections that would spoil its flatness. The aim of section 5 of this paper 
is to investigate whether the dynamics of the fibre modulus rolling down regions II and III 
gives a good description of the observed late time acceleration of Universe. 

Let us conclude this section determining the relevant scales in our set-up. As we ex- 
plained, observations constrain the dark-energy scale to be of the order pp>E ~ 10~ 120 Mp. 
In order to use the fibre modulus as a dynamical model of dark-energy, this maps into the 
following condition on the amplitude of the potential: 



F, 



poly 



(V) 



3+p 



10 



-121 



(3.12) 



This scale can be achieved by considering compactifications with very large volumes, as 
originally noted in [10], in particular for p ~ 0(1) this yields: 



(V) ~ 10 



30 



(3.13) 



in string units. We stress that, in our scenario, the VEV of the volume depends exponentially 
on the microscopic parameters of the model (see eq. (3.6)), and such large values are easy 
to obtain. So, once the minimum of the quintessence potential is set to zero, the tiny 
height of the potential, eq. (3.12), is obtained straightforwardly! Notice that values of 
the volume of the order V ~ 10 30 correspond to a string scale of the order the TeV-scale, 
M s ~ Mp/y/V ~ O (TeV), and a gravitino mass of the order the cosmological constant scale, 
7713/2 ~ Mp/V ~ O (meV), according to the discussion of section 2. Such a low- value of the 
fundamental gravity scale poses a severe problem for obtaining the original period of inflation 
since both the amount of density perturbations and the resulting reheating temperature would 
be too low. For example, the inflationary model developed in [9] requires values of the volume 
of the order V ~ 10 3 corresponding to a GUT-scale M s . The general idea that we have in 
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mind to overcome this problem is to use the volume mode as the inflaton following the work 
of [25]. Thus the volume can take smaller values of the order V ~ 10 3 in the early Universe, 
and then roll to values of the order V ~ 10 30 after the end of inflation, where the only field 
left over to roll down its minimum is the fibre modulus Tf which drives quintessence. 

A very large volume, apart from dynamically describing the smallness of the observed 
size of dark energy, has also the virtue to suppress fifth-force couplings between the quintessence 
field and SM matter, and to enlarge the very flat region II in the potential Vfibrc- I n the next 
sections, we will discuss these features in detail. 

4 Quintessence coupling to matter 

The first step in the computation of the quintessence couplings to ordinary matter is the 
moduli canonical normalisation, which we perform in detail in appendix A. The expression 
of the original moduli Tf, V, and t s in terms of the corresponding canonically normalised 
fields cj), xv an d Xs can be obtained explicitly at leading order in an inverse volume expansion. 
Considering all the moduli in the vicinity of their minimum, the canonical normalisation of 
the Standard Model modulus represented by r s , takes the form (see eq. (A. 10)): 



(Ts) 



O (V" p ) 8<p + 0(1) 5 XV + O (7v) 5 Xs , (4.1) 



where (f> is the quintessence field, while xv an d Xs are two heavy fields corresponding respec- 
tively to the volume mode and the small rigid divisor t s . The exact coefficients of (4.1) will 
be computed in appendix A at leading order in an inverse volume expansion. What is crucial, 
in the previous expression, is that the canonically normalised quintessence field <j) provides 
only a suppressed contribution to 5t s , essentially determined by the ratio of the leading and 
subleading scalar potentials of eqs (2.7) and (2.8) (recall also the general discussion of section 
2.2). Since r s controls the properties of SM physics, this implies that the coupling of the 
quintessence field to Standard Model particles is very suppressed. 

Before passing to prove this fact explicitly, let us discuss a technical issue: during the 
quintessence dynamics the two heavy fields xv and Xs sit at their minima, but the rolling 
scalar <j) can be far away from its minimum. Hence the fibre modulus has to be canonically 
normalised not just close to its minimum but for an arbitrary point in moduli space. This is a 
complicated computation whose result has already been qualitatively discussed in [26] where 
the authors used the poly-instanton effects for deriving a model of modulated reheating. 
The final upshot is that in regions far away from the minimum for large values of Tf, the 
poly-instantons get more suppressed, and so the mixing with the other moduli also gets more 
suppressed. Thus the coupling of the fibre modulus to ordinary matter gets weaker in regions 
far from the minimum. In fact, if Tf is far from its minimum then: 

e -27VT f = e -2n(r f ) e -27TT f ^ y-p^TTT/ ^ ( 42 ) 

and so the mixing term in (4.1) would be even more suppressed than V~ p . 

4.1 Coupling to gauge bosons 

We start by estimating the couplings to gauge bosons following [10]. Recalling that t s is the 
Standard Model cycle, the gauge kinetic terms are given by: 

^-gauge = —-Tr-r^r . (4.3) 
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After expanding t s around the minimum (r s = (r s ) -\- 5t s ) and performing the canonical 
normalisation of the field strength tensor, we find that the interaction term is: 

6ts_ 
IM p (t s ; 

Prom eq. (4.1) it follows that: 

0(V~ p )5(j), (4.5) 



C => 777 S ,_, G^G» V . (4.4) 



5t s 



rsj 



(rs) 
and so we find: 

£d °(m^)^ g - g " (46) 

We conclude that the coupling between the fibre modulus and the SM gauge bosons is weaker- 
than-Planckian for p ~ 0(1), in accordance with [10], and therefore poses no issues with the 
constraints from fifth-force experiments. 

4.2 Coupling to fermions 

Let us now turn our attention to the coupling between the quintessence field and the Standard 
Model fermions, which we denote as Cj. The relevant part of the SUGRA Lagrangian for 
the computation of this coupling is: 

C D KijCivfdpCi + Kjj^Hd^H + e K / 2 X i3 C i C j H, (4.7) 

where Kij and K h jj are generic functions of the moduli fields. Very little can be said 
about this theory without specifying K^ and K H j[, and so we follow the work of [27], where 
the leading modular dependence of these two functions has been estimated through simple 
scaling arguments. We simplify the set-up further by considering diagonal Kahler metrics 
and Yukawa couplings: 

C D Ki&yfdpd + K^Hd^H + e K l 2 \C i C i H. (4.8) 

Following [27] we take: 

Expanding r s and V around their VEVs as r s = (t s ) + 5t s and V = (V) + 5V we can write 
the Kahler metrics as: 

**~*^~M 1 + ^-$y) where R ^W^- (410) 

This allows us to define the canonically normalised fields Cj = y K$Ci and h = \J KqH and 
rewrite the Lagrangian as: 

5n 25V \ . „„ „„, n rr A 



CD[1 + ^T ~ ^JTTT Cii-fdyfii + dPhdph + -^rjjCiCih. (4.11) 

3(r B ) 3{V)J VKf 



L o 



Setting h = (h), expanding 1/V ~ 1/(V)(1 — 5V/(V)) and defining the fermionic mass as 
m c = ) l., /9 this becomes: 

C D clii^du + m c )ci + I —-?- — — ) c7(ry M d« + m c )ci — m c cici I . s , -\ — — | . (4.12) 
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Note that the first term is the standard Lagrangian for a free fermion of mass m c , the 
second term vanishes once we impose the equations of motion and the third term encodes 
the interaction between moduli and Standard Model fermions. Expressing the moduli in 
terms of the canonically normalised fields near the minimum, we find that the dominant 
term in the interaction with the quintessence field comes from the expression (4.1) for 5t s 
term, since to leading order 5V does not depend on the canonically normalised fibre modulus. 
This implies that the leading contribution to the dimensionless coupling between ordinary 
fermions and the quintessence field is: 

V~~o(^), (4.13) 

which is weaker-than-Planckian for p ~ O(l) like the coupling to Standard Model gauge 
bosons. 

The fact that the couplings to gauge bosons and fermions are suppressed by a scale 
higher than the Planck mass means that the current model is able to comfortably accommo- 
date the constraints from fifth-force experiments, see e.g. [3], and provide a controlled string 
theoretical description of quintessence. 

5 Dynamics of the quintessence-dark matter system 

In this section we study the dynamics of the quintessence-dark matter system by following 
the autonomous phase plane analysis originally developed in [28]. We start by analysing a 
model in which we neglect any direct coupling between the two dominant components of the 
energy density (this since their direct coupling is weaker-than-Planckian). Then we look at 
the case where there is a direct, Planck strength coupling between dark matter and dark 
energy, while at the same time dark energy couples extremely weakly to ordinary matter. In 
this way, we will outline a proposal of how to realise a coupled quintessence scenario within 
string theory. 

We consider a flat FRW Universe filled with quintessence as described by the fibre modu- 
lus and a pressureless dark matter component. The time evolution of the system is determined 
by the Raychaudhuri equation, the continuity equation for the dark matter component, the 
Klein-Gordon equation for a homogeneous scalar field and the Friedmann constraint: 

H = -^t(pdm + 4> 2 ), (5-1) 

Pdm = -3Hp DM + Q, (5.2) 

• dV Q , x 

d<p (j) 

n DM + n DE = i, (5.4) 

where O = „f.^ and Q parameterises the direct interaction between the dark matter and 
the quintessence components. Using iV = In a as the time variable and defining: 

X = AW6' ^ ^ 3MpP> (5 ' 5) 
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the equations of motion can be cast in the following form: 
x\N) = -3x 



dV/d<f> [3 y2 + 3 2 + {1 _ x 2_ y2)) + 



V 



Q 



6# 3 Mj,x' 



V'(N) = d -^\[\xy + \y{2x 2 + (1 - * 2 - y 2 )), 



H'(N) 



3H 



(2x 2 + (l-x 2 -y 2 )), 
(j)'(N) = V6x. 



(5.6) 

(5.7) 

(5.8) 
(5.9) 



In what follows we will both search analytically for instantaneous fixed points of the system 
and study the phase space trajectories numerically. 

5.1 Decoupled quintessence and dark matter 

We first analyse the case where a direct coupling between dark matter and dark energy 
is negligible, which corresponds to setting Q = 0. The system possesses two non trivial 
instantaneous fixed points, x'(N) = y'(N) = 0, whose properties are listed in table 1. These 
correspond to motion in the two distinct regions of the potential: point A corresponds 
to motion in the poly-instanton dominated part of the potential (labelled II in figure 1) 
whereas point B describes motion in the region where the loop potential dominates (labelled 
region III in figure 1). Note that the flatness of the poly-instanton potential implies that 
5 = K^oiy/^poiy ~ 0, and so the behavior of the fibre modulus approaches that of a pure 
cosmological constant. 



Point 


(Xc, Vc) 


Stability 


0,0 


^<t> 


A 

B 


(<5,Vl-<5 2 ) 

(2/a/18, y/U/18) 


stable node 
stable node 


1 
1 


-1 
-5/9 



Table 1. Fixed points for the quintessence dark- matter system. 

From the definitions of eq. (5.5), one sees that x(N) (y(N)) is the ratio of the scalar 
field's kinetic (potential) energy to the total energy. We illustrate the behaviour of the system 
in figure 2 for different sets of initial conditions. The region above the parabola corresponds 
to accelerating solutions and the dots represent the two instantaneous fixed points described 
in table 1. 

We see that if the fibre modulus starts its evolution in the part of the potential that is 
dominated by the loop effects, the system will evolve towards the fixed point B. This fixed 
point corresponds to an accelerated expansion solution, where the energy density is dominated 
by the scalar field, as preferred by current cosmological data. However the equation of state 
parameter, which is generically given by: 



<^4> 



x + y 



(5.10) 



is different from 1 for the fixed point B. Thus, given the current observational constraints, 
this fixed point is not a viable candidate for dark energy. As the fibre modulus rolls down the 
loop generated potential it eventually enters a region in field space where the loop potential is 
suppressed relative to the poly-instanton generated potential. In this region the relevant fixed 
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Figure 2. Phase plane evolution of the quintessence dark-matter system. Left: the fibre modulus 
starts in the loop dominated region of the potential; right: the fibre modulus starts in the poly- 
instanton dominated region of the potential. 

point to which the system will converge is A, which corresponds to slow roll dynamics in the 
poly-instanton dominated part of the potential. The fixed point A describes a flat Universe 
filled with a negative pressure fluid with oj = —1. Even though such trajectories display 
the same asymptotic behaviour in the future as the observed Universe, a numerical analysis 
of the system suggests they never go through a phase in which the conditions £Ide ~ 0.7 
and ujff, ~ — 1 are simultaneously verified. It then follows that trajectories that start with 
the quintessence field deep in the loop dominated region of potential do not provide a viable 
description of the late time Universe. On the other hand, if the initial conditions are such that 
the fibre modulus starts its evolution in a region of the potential where the loop effects are 
negligible, it will evolve towards the fixed point A without going through B first. As can be 
seen in figure 2, these trajectories can go through the region of phase space where the energy 
density and the equation of state parameter of the quintessence field are in agreement with 
observational data. We then conclude that there are phase space trajectories for which the 
system evolves to a state where it provides a viable description of an accelerated expanding 
Universe, namely the ones corresponding to initial field values in the poly-instanton plateau. 
It is worth noting that from the discussion in section 3 the width of the poly-instanton 
plateau is determined by the ratio V^, i y /Ui oop ~ V 1_p which can easily be > 1 given the 
large volumes required for quintessence phenomenology, V ~ 10 30 . There is therefore an 
appreciable range of initial conditions for the fibre modulus that yield a viable quintessential 
behaviour. 

When looking at the late time composition of the Universe one is confronted with the 
coincidence problem: why is the dark energy density comparable to the dark matter density 
today? Dynamical models of dark energy often address this issue by allowing for scaling 
solutions: fixed points where the ratio of dark energy to dark matter density is constant. 
The existence of these scaling solutions is determined by the steepness of the dark-energy 
potential (V /V) and by the equation of state for the matter component. In a Universe with 
pressureless cold dark matter and the fibre modulus playing the role of quintessence, there 
are no such fixed points as both regions of the quintessence potential are too shallow to 
support them. Therefore we are led to the conclusion that even though our model provides 
a realistic description of the current state of the Universe, it does not allow for a dynamical 
explanation of the coincidence problem: we just happen to be around at a time when Qde 
and &dm are of the same order and the Universe will eventually evolve to a state where it is 
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filled with quintessence. 

In these considerations about the distant future evolution of the system we have so far 
ignored the presence of the minimum of the potential, region I in fig. 1, to the left of the 
poly-instanton plateau. Its presence implies that in the very distant future the fibre modulus 
will eventually fall into its minimum, after which will follow a phase of oscillatory dynamics, 
in a process similar to post-inflationary reheating, albeit at a much lower energy scales. Since, 
ultimately, acceleration ends with the scalar sitting at its Minkowski minimum, we avoid the 
problems raised in [29, 30] with respect of embedding asymptotic de Sitter space in string 
theory. 

5.2 Coupled quintessence and dark matter 

We now analyse how a coupled quintessence-dark matter system affects the dynamics of the 
Universe [31]. Possible Planck suppressed couplings between dark energy and dark matter 
are allowed in our set-up if dark matter is constituted by bulk KK modes or moduli fields. 

Let us start discussing possible dark matter candidates: phenomenological constraints 
on dark matter from collider physics, astrophysics and cosmology in this class of scenarios 
have been discussed in [10] . Here we briefly summarise the main conclusions. The requirement 
of avoiding BBN constraints, associated with the fact that the SM brane cools too quickly 
due to evaporation of energy into the bulk, sets a severe bound on the reheating temperature: 
2rh ;$ 100 MeV. If this constraint is satisfied, there are no issues with BBN, but the relic 
KK modes, if stable, can still overdose the Universe. In order to avoid this, the previous 
bound becomes more severe, bringing the reheating temperature at most around a few MeV, 
just above the BBN temperature. On the other hand, if the KK modes are not stable and 
decay into photons, then one can have problems with possible distortions of the diffuse MeV 
7-ray background. 

Notice that these bounds on the reheating temperature imply that we cannot have a 
thermal bath up to the weak scale, and so the standard thermal WIMP picture for dark 
matter does not work. In addition, we do not have neutralini in the 4D low-energy spectrum 
since supersymmetry is badly broken on the SM brane. In our case we have two kinds of 
KK modes: 6D KK modes and 10D KK modes. The 10D KK modes decay whereas the 6D 
KK modes are so long lived that can be considered as almost stable with respect to the age 
of the Universe. Hence given that the 6D KK modes are almost stable, we need to have a 
reheating temperature just above the MeV scale. The advantage is that the 10D KK modes 
decay predominantly to 6D KK modes, and so their decay to photons is suppressed. Thus 
we have no problems associated with 7-ray background distortion. 

Higher 6D KK states can be produced from the thermal bath but they quickly decay 
to lower 6D KK states, which are then long-lived. There are also two light moduli which are 
stable with respect to the age of the Universe. One is the volume mode which we can suppose 
to be the inflaton, so that inflation takes place when the gravity scale is high, whereas at the 
end of inflation, the string scale is around the TeV. This implies that the volume mode is 
very light with a mass of the order meV. Hence this modulus can be produced in the thermal 
bath (or by simple scalar oscillation) and then suffers from the cosmological moduli problem 
(CMP). The other light modulus is the fibre divisor which drives quintessence. 

Therefore, we conclude that in our scenario dark matter can be realised in terms of a 
mix of bulk 6D KK modes and the volume modulus, once the CMP associated with the latter 
is solved (for example via a second low-energy inflationary period as in thermal inflation). 
The quintessence field would then indeed couple to dark matter, for example 6D KK modes, 
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with Planck suppressed trilinear metric couplings of the form: 

i«W, (5.H) 

where <j> is the quintessence field and ip the 6D KK mode. 

After this qualitative discussion to motivate the coupling of dark energy to dark matter 
in our set-up, we can write the dark matter energy density as a function of the fibre modulus 

as: 

PDM = Po/O), (5.12) 

where po obeys the normal conservation equation: 

po = -3Hp , (5.13) 

and f(4>) parameterises the interaction between dark matter and the quintessence field. The 
conservation equation for dark matter is then given by: 

PDM = -3Hp DM + Q , .fjs pDM- ( 5 -14) 

By comparing eqs. (5.2) and (5.14) we define Q = —jMlTJiPDM- I n order to look for 

fixed points analytically, one must be able to write tt^qj in terms of x(N),y(N) and other 
dimensionless parameters. One simple case might be: 

/(</,) = e a ^ Mp , (5.15) 

since this form for the coupling function implies that: 

wm) PDM = ^ PDM ' ( } 

allowing for an analytical study of the phase plane dynamics of the system. This is exactly 
what one expects to find if the dark matter mass is a function of the fibre modulus r/, 

m DM ~ T f „ e 2mMV3 ( 5 _ 17) 

where in the last step we have written r/ in terms of the canonically normalised field (f>. 
It then follows that Q = —a\^6xHp£>M and eq. (5.6) becomes: 

X > {N) = _ 3x _ ^^y|y 2 + \x(2 X * + (1 - x 2 - y 2 )) - a^l - x 2 - y 2 ). (5.18) 

The phase space for models with this form of the direct coupling Q has been studied 
in detail in [32], where a phenomenologically interesting fixed point located at (x c ,y c ) = 
(^f , V9 + 6 fc afe ) was found. Note that we have defined A = -M P ^ and b = A + a. Following 
[32] we find that this fixed point exists if: 

3 3 

b 2 > - and - - < ab < b 2 - 3 (5.19) 
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and is stable if: 
and 



,2 3 
b 2 > - 



- < ab < -B sgn ( b ) 



or 



B_ sgn{b) <ab<b 



where: 



B± 



2b[±(b 2 - 3/2) 3 / 2 - b(b 2 - 39/8)] 



(5.20) 



(5.21) 



46 2 + 3/4 

As in the decoupled case, we proceed by considering two different regimes, depending 
on which term in the potential dominates. For each case we shall then investigate whether 
the above conditions are met. In the regime in which the poly-instanton effects dominate 
over the loop potential (region II of figure 1), A ~ to a very good accuracy. One sees that it 
is impossible to satisfy the existence conditions given by eq. (5.19). We then conclude that 
the poly-instanton potential is unable to support an interacting tracking solution. However 
if one looks at the regime where the loop generated potential dominates (region III of figure 
1), one has A = — 2/y3 and one sees that there are values of a for which a stable accelerating 
tracking solution exists. This region however is not very wide due to the requirement of 
stability. In fact we find that it is roughly —-^-m — 0.55 < a < — 7775- We must note that 
this is not so much a consequence of the shape of the potential as of the explicit form of Q 
we are assuming. Furthermore, we see that in this narrow stability region 0.93 < Qde < lj 
where the lower bound corresponds to a = — -^= — 0.55 and the upper bound to a = — -A=. 

2v3 2V3 

It is worth noting that for a < the dark matter mass goes like tudm oc I/tJ 1 . 

The behavior described above is illustrated in fig. 3, where it is clear that the presence 
of the interaction dramatically changes the trajectories in phase space. One also sees that 
the coupled system (magenta) converges to a fixed point which is different from the one of 
the decoupled case (blue lines). 
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Figure 3. Phase plane evolution of the quintessence dark-matter system with interaction and -F po iy 3> 

-^ loop 

Once one considers the full potential, the fixed point described above is transient since 
as the field rolls down the loop generated potential it eventually enters the poly-instanton 
dominated region. So the phase-space trajectories start by converging to the fixed point B 
before evolving to the cosmological constant fixed point A. 

As in the non interacting case, the fibre modulus is able to provide suitable description 
of dark energy even though it does not address the coincidence problem. The introduction of 
the non minimal interaction term extends the existence of viable trajectories for cases when 
the field starts its evolution deep in the loop dominated part of the potential. 
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Figure 4. Phase plane evolution of the quintessence dark-matter system with interaction and i^ po iy ~ 
10 -Fi 00 p 

6 Conclusions 

In this paper, we presented a new quintessence model embedded in string theory where the 
small mass of the rolling scalar field is technically natural and its coupling to Standard Model 
degrees of freedom is naturally weaker-than-Planckian, thus avoiding fifth-force constraints. 
We focused on a type IIB LARGE Volume Scenario with a very anisotropic shape of 
the extra dimensions since two of them are stabilised at values exponentially larger than the 
other four [10]. This leads to a string derivation of SLED scenarios [11] which have already 
turned out to be very promising for obtaining models of quintessence [12]. In these scenarios, 
if TeV-scale gravity provides the solution to the hierarchy problem, 4 the gravitino mass m 3 / 2 
turns out to be of the observed order of magnitude of the cosmological constant: m 3 / 2 ~ A. 
Therefore, the current value of dark energy A can be naturally reproduced if the quintessence 



field develops a potential at order V 



~ m 



3/2' 



This mechanism, which ties the solution of the hierarchy problem to the one of the 
cosmological constant problem, relies on the existence of a compensation between the SM 
brane tension and its back-reaction. A possibility to achieve this could be to directly apply 
to our set-up the ideas developed in 6D SLED constructions. This, in turn, depends on the 
details of the 6D effective supergravity theory for our anisotropic compactifications; on the 
other hand, this intermediate theory has not yet been explicitly identified in [10]. We leave 
a proper understanding of these interesting open issues for study in the near future. 

In this paper, we showed that the fibre divisor T/ indeed develops a potential at order 



V 



~ m 



3/2 



via tiny poly-instanton corrections to the superpotential and closed string loop 



contributions to the Kahler potential. Furthermore, this field is radiatively stable since the 
bulk is nearly super symmetric. 

We also showed that our quintessence model can evade the severe bounds coming from 
fifth-force experiments even if we are dealing with a scalar. The reason is the fact that the 
rigid four-cycle supporting the Standard Model does not intersect the fibre divisor, and so 
the quintessence scalar couples to ordinary matter only via its mixing with the volume mode 
which, in turn, induces its mixing with the SM cycle. However due to the fact that the 
fibre modulus is a flat direction at leading order which is lifted only at subleading order, the 



4 Notice that all the scales of our model are set by the value of the overall volume V which has an exponential 
dependence on the underlying parameters, and so it can be easily set to arbitrarily large values without tunings. 
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mixing between tj and V is suppressed, resulting in a weaker-than-Planckian coupling of the 
quintessence field to SM particles. 

On the other hand, if dark matter is realised in terms of quasi-stable Kaluza-Klein states, 
direct couplings between dark energy and dark matter are allowed, leading to a scenario 
of coupled quintessence within string theory. We studied in detail the dynamics of the 
quintessence scalar in our set-ups, investigating the nature of fixed points and the late time 
evolution of dark energy, showing that the main features of our scenario make it compatible 
with observations. 

To conclude, our model explicitly shows that concrete string theory constructions are 
able to lead to scalar theories of quintessence that, besides generating naturally flat scalar 
potentials, also provide couplings to Standard Model matter that are weak enough to satisfy 
present constraints on fifth forces. 
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A Canonical normalisation 

In this section we compute the canonical normalisation for the fields involved, using the 
well-known techniques developed in [27, 33-35]. The Lagrangian for the quadratic field 
fluctuations around their minima, T{ = (tj) + <5tj, takes the form: 

c = Kijdn&ndPSn - (V) - ^Vij&n&rj , (A.i) 

where Kij is the Kahler matrix, and Vij the matrix of second derivatives of the potential. 
The leading terms in the Kahler matrix Kij look like: 



/ 3 _J 3y/^ \ 

8tJ 4Vr f 8Vr f 



K ~ 



i i 



4Vt/ 8V17 

1 _vi 3 I 



(A.2) 



where we used (tj, V,t s ) as coordinates of the Kahler moduli space. We write the original 
moduli in terms of the canonically normalised fields around their minima as: 

Sr f \ fSd>\ 

SV \ = C \8xv , (A.3) 

St s J \$Xs J 

where C is the matrix that diagonalises both the kinetic and the mass terms. It can be 
written as: 

C = ( v {1) v {2) v {3) ) , (A.4) 
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i v-h 



where vu\ are the eigenvectors of the mass matrix Mj/% = ^ K- Vjk, normalised such that 



Cj k Kkm,C m j -- Sij. The mass matrix reads explicitly: 



Af, 



lJ ~ 2 



4t / (/ 8 +2/ 6 t / +2/ 7 t s )V 4(/ 3 +2/ 4 )r /v ^7 4(/ 4 +2/ 5 )r / \ 

V3+p(3r s 3/ 2 -2V) (3r s 3/2_2V)V3 3r s 2 V 2 -2 v ^V 3 

2V(2(/ 6 r / +/ 7 r ;i )V+3/8(-r s 3 / 2 +V)) 6/ 3 r s 2 -2(3/3+2/ 4 ) v ^V 6/ 4 (r s 3 / 2 -V)-4/ 5 V 

V 3 +p(3t s 3/ 2 -2V) (3t s 3 / 2 -2V)V 3 3r s 2 V 2 -2 v ^7V a 

4(3(/8+2/ 6 r / )T 3 V+4/ 7v ^V 2 ) 4(3/ 3 r s 3 / 2 +4/ 4 V) 4(3/ 4 r s 3 / 2 +4/ 5 V) 

\ V 3 +p(9t s 3 / 2 -6V) 3V 3 (-3r s 3 / 2 +2V) 3V 2 (3r s 5 / 2 -2r s V) / 

(A.5) 

The quantities /$ are independent from the volume at leading order in an inverse volume 
expansion. Their expressions are quite long, and we relegate them in appendix B. In writing 
the previous mass matrix, we added to the leading part of the potential, eq. (3.4), the 
subleading contribution due to polyinstanton corrections, eq. (3.5). 

Notice that, as expected, the first column of the mass matrix is suppressed by a factor 
1/V 3+P , since it is due to the subleading poly-instanton potential. In other words, the first 
column of the matrix (A.5) would be zero in the absence of the subleading poly-instanton 
contribution. 

It is straightforward, but long, to extract the eigenvalues and the eigenvectors of the 
mass matrix of eq. (A.5), working at leading order in an expansion in inverse powers of the 
volume. The expressions of the eigenvectors turn out to be particularly complicated. It is 
convenient to present them at leading order in an expansion in the quantity r s , that we can 
assume to be small with respect to the other quantities. We find: 



C 



( *n A /2 r . \ 

y/S V 3 f 



ci 



V6i 



gTsVVv 



(A.6) 



\ W Vur s ^ / 

The two last entries in the first column are suppressed by powers of the volume, and are due 
to the effect of the poly-instanton corrections. The quantities cq and c\ are independent from 
the volume (at leading order in a inverse volume expansion) and read: 

co = -^3-l^-2aAC 1 , + bB { c(-C 1 + C 2 ) + 2C 1 ,)) (^^) ^ (A.7) 

4r| 3-l + P(2aAC 1 n + bB(c(Ci - C 2 ) - 2C l7 r)) (^^Y 
Cl = (^^ ^ L (A-8) 

The lightest eigenvalue of the mass matrix scales as m? oc 1/V 3+P , and provides the effective 
mass for the quintessence field in the poly-instanton dominated region of the quintessence 
potential (region II of figure 1). 

In this region, the volume and t s are stabilised at their minima since their mass is very 
large. Prom eqs (A. 3) and (A.6), one obtains that the displacement of the light field r/ from 
its minimum can be expressed in terms of its canonical counterpart as: 

2 _2_ . 

St/ cz —j= Tf 5<p =>■ Tf ~ e v^ , (A. 9) 

v3 

justifying the canonical normalisation we used in eq. (3.8). The diagonalising matrix (A.6), 
when compared with (A. 3), also tells us that, within the approximations we are using, we 
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can write the modulus r s in terms of canonically normalised fields as: 



St. 



c\ 



(t s ) " VP 
reproducing the scalings of eq. (4.1 



T==&<t> + yfc8xv + 



2W 



V3(r s )3/4 



$Xs 



(A.10) 



B Explicit form of the coefficients of the mass matrix 

Here we collect the complete expressions for the quantities fi appearing in the mass matrix 
Mij, eq. (A. 5) (we use the same conventions introduced in our section 3): 



h 

h 
h 



((a - b)bBr s + Z- ar s Z) 
(4(a - b)bBr s + Z - Aar s Z) 
3VWo (b 2 Br 3 -Z + aT s (-bB + Z)) 
Z (Ab 2 Br s - Z + 4aT s (-bB + Z)) 

| (-36(-l + h)fiT s W Q 2 - 9f 2 ^F s W Z + 2f 2 2 Z 2 ) 



U = 3/2 (-/2^(4(a - b)bBr s + Z- 4ar s Z) + 6^W ((a - b)bBr s + Z- ar s Z)) 
/5 = 0/2 



6T s 3/2 W (b(-a + b)B(-2 + (a + b)r s ) + o(-2 + ar s )Z) + f 2 ( 8(0 - b) 2 b 2 B 2 T 2 
+8b(-a + b)BT s (-l + 3clt s + 6r s )Z + (-1 + 8ar s (-l + 2ar s ))Z 2 ^ 



h 



16tt/. 



1+p 



3p 



4//+P r 



36B(-c(Ci - C 2 )(-l + ttt/) + (-ad + bC 2 )iTT s )W 
+vr (4B(oCi + 6((-l + b)Ci - bC 2 ))f 2 ^r~ s + 3Ci(-2 + 2^r/ + or s )W ) Z - 4aCi/ 2 vr^Z 2 
f 2 



3e./rT 



365^ {-2aCm - 2b 2 C 2 iTT s + &(c(d - C 2 )(-l + 27TT/) + 2tt(C 2 + aCrr,))) W 

-2tt (2B(ad + 6((-l + 6)Ci - 6C 2 ))/ 2 (-l + 46r s ) + 3Ci^(6(-l + 27TT/) + a(-l + br s ))W ) Z 
+4aC 1 f 2 ir(-l + 4bT s )Z 2 ~ ] 



fs 



i+P 



3e 



36 J B(c(Ci - C 2 )(-l + 27TT/) + 2(oCi - bC 2 )itT s )W Q 



-2vr (2B(aCi + &((-! + b)C x - bC 2 ))f 2 ^rS + 3C*i(-l + 27TT/ + ar s )^ ) Z + 4aCi/ 2 7rV^Z 5 
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